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1. INTRODUCTION 

As one of promising approach implemented in 5G systems, non-orthogonal multiple access (NOMA) 
was introduced |I]. To provide massive connections, NOMA can rely on this first main benefit to further serve 
services with high spectral efficiency and low latency. Those advances are urgent requirements to design new 
generation of 5G and beyond wireless systems [2]. In NOMA, the multiple users can be shared same frequency 
but different power levels are assigned to each user effectively. The signal detection technique is required at 
receiver to extract information exactly with low error. How NOMA treats far users and near users to assign 
power levels. Fortunately, by detecting the channel gains of different channels, suitable power coefficients 
are assigned to users reasonably [3)—[5]. As interesting application of NOMA techniques, half-duplex relay 
stations (RSs) have been studied in order to increase spatial diversity [6J—[13]. The benefit of relay can be 
reported in [8], [T0], and [OT], since single relay is placed between transmitters and receivers. Nakagami-m 
fading channels and Rayleigh fading channels are popular channel models deployed in the NOMA 
system relying on a single amplify-and-forward (AF) relay, which outperform the orthogonal multiple access 
(OMA) in terms of two system performance metrics (outage probability and throughput). Kader et al. in 
developed a network containing two sources, two destinations, and a relay to form the cooperative NOMA 
with a half-duplex decode-and-forward (DF). They examined perfect and imperfect successive interference 
cancellation (SIC) when they evaluated ergodic sum capacity. 

As simpler approach, Liang et al. and Xu et al. in [I4], studied half-duplex (HD) relay-based 
NOMA systems. Due to the requirement of additional time resources, HD NOMA just provides low spec- 
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tral efficiency. Different from HD NOMA, full-duplex (FD) relay-based NOMA deploys the same frequency 
channel to permit the relay to simultaneously receive and transmit signals and reduce such loss [16]-[19]. The 
operation of FD benefits from advances of antenna isolation and cancellation of analog self-interference (SI) 
at the FD relay. FD device-to-device assisted coop erative NOMA system was investigated in which the 
near user needs the FD relay to support transmission to the far user. Deng et al. adopted Rician fading 
channels for FD NOMA system by evaluating formulas of the outage probability and the ergodic rate under 
imperfect conditions such as imperfect SIC and residual hardware impairments at transceivers considered 
the impact of imperfect SIC and residual inter-relay interference on a DF relaying based NOMA. The authors 
developed for the considered framework over generalized Nakagami-m fading channels by evaluating outage 
probability (OP), asymptotic OP, and ergodic rate. The work in studied downlink NOMA short-packet 
communication systems the average block error rate (BLER) by using stochastic geometry and Nakagami-m 
fading channels. A few work consider FD at relay for NOMA, for example [24], which motives us to study 
difference among two destinations under the impact of hardware impairments. 


2. SYSTEM MODEL 

A dual-hop NOMA transmission with the help of a FD relay (R) is studied, shown in Figure.}1} The 
system model could be examined in the case of a base station (B) serves a dedicated group of two NOMA users. 
In particular, we design a FD relay (R) which is intermediate device while two NOMA users including D; and 
Dg. Those users are classified based on channel gains to determine the near and the far users. The FD relay is 
equipped two antennas to transmit and receive signals simultaneously. To provide general channel model, all 
the channels are assumed as Nakagami-m channels. We treat power coefficients £1, €2 to help the base station 
B serving dedicated group of users and satisfying strict constraints ¢, + €2 = land £1 > €9. 
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Figure 1. Considering hardware impairment aware FD NOMA system 


The transmit signal processed at the base station B is /é,P p21 + VéoP B22, in which Pg is the total 
transmitted power of B; x1, £2 are denoted as signals of D1, Do. 

The transmit signal from the base station is then processed at the FD relay. It is worth noting that FD 
relay produces xy, as the loop self-interference which is expected to eliminate. In particular, we can compute 
the received signal at the relay as: 


yr = go (yB + ns) + gré (v Prig, +m) + No; (1) 


where go, gf are the channel coefficients of B + R and R —> R links. € is denoted for HD/FD modes, i.e. € = 0 
and € = 1 are known as HD and FD modes respectively. ng ~ T (0, KeRPBIg9o| jea m~r (0, Ke Prigs| g: 


n; ~T (0, Kp, Prigil *) , (i = 1,2) represents noise distortion, no ~ I (0, No) denotes the additive white 


Gaussian noise. KBR, Kg, and Kp,, (i = 1, 2) are the levels of residual hardware impairments. 
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The second hop signal processing is conducted based on channels g; which is refereed to links R —> 
D;. The NOMA power allocation need be adjusted for the second hop transmission, i.e. 41, p2 are allocated to 
two signals of two users, those factors are satisfied 41 + u2 = 1 and py > u2. 

The received signal at Dı, Də are given as: 


YD; = Ji (v PRT + y H2PRT2 + no: ) + no, (2) 


The two signals x1, £2 need be processed at the FD relay based on signal to interference plus noise 
ratio (SINR) which can be computed by: 


2 
€1p8|90| 
Ya1,R = 2 2 Hi (3) 
(€2 + KBR) PBlgol + (1+ KBR) PREP lgsl +1 
and 
ce e2pBlgol” (4) 
£2, R — 2 2 , 
pBKRRIGo| + (1+ Kr) PRE If] + e1pBhR +1 
where pg = 5B, PR = ie are the transmit SNR at B and R. ha ~T (0, w\g0l”) caused by imperfect SIC 


(ipSIC) and w € [0, 1). 
After signals transmitted at the second hop transmission, the destination Dz need to know SINR as 
below. In particular, Dz detects signal xı as shown in: 


2 
H1 Pr|92| 
Yas,Da = PR ee — e, (5) 
L2prig2| + prr, lg2l +1 
The first user D; wants to detect signal x1, x2 respectively as shown in: 
mprlgl? 
Yxı,Dı = 2 2 2 ’ (6) 
Hoprigil + prkp,lgail +1 
and 
uprlo? 
Ya2,D1 = 2 2 bj (7) 
HıprhrRD, + prKp, lgl +1 
where hrp, ~T (0, wlg1?). 
The PDF of the Nakagami-m channel gain gx (k = 0, 1,2, f) can be expressed as: 
fog? (2) = ee Pe (8) 
2 (x) = ———,€ lr, 
lgx l T (Mg, ) N 
where bg, = Aarf is the mean value of gẹ denoted by | g|? ~T (trisi dop T The CDF can be written by: 
Ik Mygg 
Mg), —1 
1 x — gr 
F 2(2) =1- pp (mo) =1-e Pox : (9) 
lgx r (Mgr) s Bor 2 npr. 


3. PERFORMANCE ANALYSIS 
3.1. Outage probability of Dı 
To evaluate system performance, the OP need be compute, the OP of D; is defined as [19], [20]. 


min aes > YP, Yaar 2 AP), ) 


OPp, =Pr (min Za, < 7t)) + Pr : 
ve ( er TA aR i )) min (Yzı,Dı < YAP, Yez, D1 < vin) 


(10) 
= 01 +02 
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where the threshold SNRs are yi? = 21 — 1, 8h = 2®2 — 1. 
Replacing the formulas (3), into (LO). we can calculate 9; as shown in: 


x path 2 th path 2 th 
6 Spr gol? < EE E jg Ay) Pr op < SHA YW jg A 
pı pı p2 p2 (11) 
= m1 +T, 


A A A A 
where yı = €1pp — (E2 + Kip) PBIS, G2 = ezpp — (E1w + KÊR) PB}, ® = (1+ Kh) PRE, Ar = 


aate and we can be calculated T1, T2 as: 
2 


th 
Ti 2 Pr lg | Y2 lorl ae Jigar < Ai 
Yl 
Ai Mg,—1 
af at E 90 path th ” 
= x e Bay dz — 5 Ya T Ea avis ayant tedy (12) 
I (mgr) bor” pı pı 
gf of n=0 69 


2 T 
i A ?18 9 A py 1 e 2890 A oy 1 
In which = e = 2 a = 1 a 
n which, œi nian (ma) bo » Q2 P1899 Bap’ 5 nlar (ma Ba > A6 P2899 Bag 
and T> is computed by. 
T Ap, lgo]? < 7i "gs? 2> A, 
p2 Tas 
oo 1 
a Baz TE n ai ” k „ktmMmg,—1 eT oer 
= | st da LX S(p (Ë) asatattmytersorde a3 
A L (Mgs) Bas n=0 k=0 
1 
ma, Magy — n y n -k-mg; ‘ 
= Bg, T mass - > 5 A5 AG ® T ((k+mg,),a6A1). 
n=0 k=0 p2 
Similarly 01, 02 can be calculated as shown in: 
A : : 
Oo S Pr (min (aa > 72"; Yoa,R = 71") Pr (min (Ye1,D. < 73°) Ira, < 71") (14) 
= T3 X T4. 


We can calculate 73 and 74 as shown in: 


A . th yte 
T3 21 — Pr (lgo]? < min ( Asa gy k 47 i, Aer lor lg € )) 
= 1- yı- Y2 


(15) 


A A A 
where Ay = €1pp — (£2 + KBR) PB, As = e2pB — peKBRII’ — E1pBwyit, Aa = (1+ rir) PRE’, 
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13 
as A (Ey e 42899 A‘ o ê _ A271 —Asys" is A Aas 1 
Az nlp T (mg, ) Bae? Í AsAayg?—AoAar”? A2Bgq Bas” 


A 2 Asy 2 yi 2 
v Pe (lo < SE loy + B, lg? <0 
2 2 


O mgo—1 n 


© 
Mg, —1 
P If eee n 2 _ 
= | ——— ne °F dx — J ) agr tas leat gy 
I (m ) If kı 
3 gs) Pos 6 n=0 ky 
O Mg —1 
ei (mo), E) g0 n 


= m —ki-mg, 
T (mg;) 2 mt ky Jai azy ((kı + Mg;) ,048) , 


=0 


A 9 Aate 9 yte > 
Wo 2Pr (lw < A, lor t Ap las 20 
e = j i 
_ T Gar 2) _ y y ( n ) (£) Abe Ket 
I (mg) n=0 ky kı As mpg. (Mgs) Bas? 


xT (r Lm A47 O _ 2 ) Gece o1 jo 
an? A3 bgo Bags A3 bgo Bos ` 


Then, 74 can be expressed as: 


A . pıpr|g |? th u2prlg1|? th 
=Prin < < 
TA r ( m ( w2prigil tere, |g? +1 V2 mprolg| teres, l N 


= Pr (Io? < min (As, As)) 


-1 
_min(4s.46) "ZLU (min (As, Ag))” 


=l-e Bay 
Qn 
n=0 n! 
S A 4th A ath 
with As = = E 2 „ths A6 = = LEL -2 nth: 
H1PR—H2PRY2 -PRED V2 H2PR—HMIPRWY, -PRÉD V1 


(16) 


(17) 


(18) 


We have applied the formulas [25], (1.111), [25], (6.381.1), and [25], (3.381.3) in the calculation steps 


above. 


3.2. Outage probability of Də 
The OP of D2 can be written as: 


OPp, = Pr ( Yzı,R < Y3, Yar sls < WP, ear < yo ) 


= Pr (Yar < YP, Yaz, R < vin) x Pr (Ye1,De < ait) 
= 0, x 03, 


where 0; was calculated in the previous section and afer substituting (3) into (19). 03, we obtain: 


th 
A 2 Y 
83 =Pr | |g2|° < th 2 ~th 
H1PR — H2PRY2 PRED, V2 


yir ) a 
th m 1 2 
72 g1 th 2 th 
= = K 
=J—e "1PRÊI1-#2PRÊ91 15? — PRED, Bg) 75" 5 Bape FRPR PRRD 12 
n! Br 
n=0 gı 
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4. SIMULATION RESULTS 


In this section, we assume that the levels of RHIs K = KBR = Kg $ KD, KD, the mean values 
of channel power gains Ag, = Ags, Agi Ags» the target rates of Dı, De are respectively Ri, Ro, w = 0.01 
and power allocation coefficients €; = H1, €2 = u2. The better quality of channels (higher m) leads to 
improvement of OP performance for two users, shown in Figure 2. In addition, in Figure 3, higher requirement 
of data rate R4, Rə results in worse OP performance. The reason is that in (TO), OP depends on the target rates. 


We then see the impact of level of self-interference channel at the relay on OP in Figure 4 performance. 
Ag; = 0.1 is reported as the best case for two users. The difference among two users is decided by different 
power allocation factor assigned. The impact of hardware impairment can be observed in Figure 5. Less impact 
of hardware impairment «x = 0.001 is the best OP performance. 
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Figure 2. Outage probability versus transmit SNR (p) Figure 3. Outage probability versus transmit signal to 


with different m with R, = 2 (bits/s/Hz), Rə = 1 noise ratio (SNR with different target rates, 
1 2 
(bits/s/Hz), €&1 = py = 0.6, k = 0.05, Ago = Àg = 8, E1 = fy = 0.7, k = 0.05, Ago = Ags = 5, Ag, = 1, 
Ag = 1, Ags = 0.5 Ags = 0.5,5m=2 
10° 
0.7 4 
0.6 3 
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Figure 4. Outage probability versus transmit SNR (p) Figure 5. Outage probability versus transmit SNR (p) 


with different of Ag re Rı = 0.5 (bits/s/Hz), Rə = 0.5 with different of K, Ry = 2 (bits/s/Hz), Rə = 1 
(bits/s/Hz), Ago = Ago = 5, Ag, = 1, Mm =3, (bits/s/Hz), €1 = uı = 0.66, Ago = Ago = 8, Ag, = 1, 
E1 = u41 = 0.6, k = 0.01 Àg; = 0.5,m =2 
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5. CONCLUSION 

In this article, a downlink FD NOMA system was studied under the impact of hardware impairment. 
To illustrate advantage of NOMA scheme, the closed-form expressions of outage probability were provided. 
Numerical results were presented to corroborate the theoretical analysis, demonstrating that the quality of 
channel, level of hardware noise yield significant performance gains over Nakagami-m fading. Moreover, all 
the results showed that the system performance is limited by the target rates. NOMA with more users can be 
addressed in the future work. 
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